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Abstract 

The  Klein-Nishina  formulae  for  the  differential  and  integral  cross  sections  for  Compton  (or 
incoherent)  scattering  of  x-ray  photons  from  atomic  electrons  assume  that  the  electrons  are  free 
(unbound)  and  at  rest,  which  is  a good  approximation  for  photons  of  the  order  of  1 MeV  or  higher, 
particularly  for  low-Z  target  materials.  The  probability  for  dislodging  bound  electrons  can  be  taken 
into  account  at  least  approximately  by  including  an  incoherent  scattering  function  S(x,Z)  in  the 
differential  Klein-Nishina  formula.  Currently  available  systematic  theoretical  tabulations  of  S(x,Z), 
where  x is  a momentum  transfer  variable  dependent  on  the  photon  deflection  angle  and  the  photon 
incident  energy  do  not  use  S-matrix  theory  with  its  claimed  higher  accuracy.  The  currently  most 
widely-used  S(x,Z)  tables  are  non-relativistic.  Recently  (1998)  Kahane  has  used  the  Ribberfors- 
Berggren  relativistic  treatment  to  compute  systematic  tables  of  S(x,Z).  These  new  values  of  S(x,Z) 
have  not  yet  been  exploited  for  computing  relativistic  bound-electron  incoherent  scattering  integrated 
cross  sections  a^,  nor  for  triplet  (pair  production  in  the  field  of  the  atomic  electrons)  cross  sections 
Ke  which  also  require  S(x,Z)  for  their  computation.  The  motion  of  the  atomic  electrons  around  the 
atomic  nucleus  gives  rise  to  a Doppler  broadening  of  the  apparent  energy  of  the  incident  photon, 
resulting  in  a corresponding  broadening  of  the  Compton  "modified  line"  for  a given  deflection  angle 
of  the  outgoing  scattered  photon.  The  shape  of  this  broadened  line  is  called  the  "Compton  profile." 
One  task  of  this  report  is  to  examine  the  usefulness  and  advantages,  if  any,  of  S-matrix  theory  in 
providing  significantly  more  accurate  values  of  photon  interaction  cross  sections  and  attenuation 
coefficients.  The  other  task  of  this  report  is  to  examine  the  available  Compton  profile  literature  and 
to  explore  what,  if  any,  effect  our  knowledge  of  this  line-broadening  has  on  theoretical  computations 
of  photon  incoherent  scattering  cross  sections  and  total  mass  attenuation  coefficients. 
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1.  Introduction. 


S-Matrix  Incoherent  Scattering  Theory: 

The  availability  of  S-matrix  theory  for  computing  photon-atom  interaction  cross  sections  was  briefly 
touched  on  in  a previous  report  in  this  series  (Hubbell  1997).  This  report  will  now  endeavor  to 
estimate  the  magnitude  of  the  anticipated  impact  and  possible  adjustments  in  values  of  the  mass 
attenuation  coefficients  ji/p  from  those  in  the  above  compilations,  in  changing  the  theoretical  bases  for 
systematically  computing  the  incoherent  scattering  cross  section  from  those  of  Hubbell  et  al 
(1975)  based  on  the  work  of  Brown  (1970a, b,  1972,  1974),  of  Cromer  and  Mann  (1967a, b)  and  of 
Cromer  (1969,  1971)  to  the  S-matrix  approach.  The  S-matrix  approach  is  outlined,  for  example,  in  a 
series  of  papers  by  Bergstrom  et  al.  (1992,  1993a,b,  1995,  1997)  and  by  Pratt,  Kissel  and  Bergstrom 
(1994).  In  the  survey  in  Section  2,  and  in  the  conclusions  in  Section  4,  the  feasibility  and 
implications  of  S-matrix  incoherent-scattering  theory  for  computations  of  the  Compton  and  triplet 
(Hubbell  et  al  1980)  contributions  to  the  energy-absorption  coefficient  \lJp  (Seltzer  1993)  is  also 
explored. 

Available  Data  on  Compton  Profiles: 

Since  the  impact  of  Compton  profile  data  on  the  computation  of  x-ray  cross  section  and  attenuation 
coefficients  seems  not  to  have  been  explored  previously  in  the  literature,  the  major  part  of  this  report, 
Section  3,  will  focus  on  this  second  topic. 

In  the  early  experimental  results  by  Gray  (1920),  and  by  Compton  (1922,  1923)  on  the  x-ray 
scattering  effect  bearing  his  name,  and  in  the  theoretical  work  by  Debye  (1923),  by  Dirac  (1926),  by 
Klein  and  Nishina  (1929),  by  Compton  (1935),  and  later  by  Platzman  and  Tzoar  (1977),  little 
attention  was  given  to  the  apparent  broadening  in  energy  of  the  deflected  photon,  due  to  the  poor 
resolution  of  the  early  spectrometers.  Some  early  theory  was  provided  by  Jauncey  (1925a,b,  1934), 
after  which  some  definitive  measurements  of  this  broadening,  or  "Compton  profile",  were  made  by 
Du  Mond  (1929,  1930,  1933)  and  Du  Mond  and  Kirkpatrick  (1931,  1937)  using  an  x-ray  tube  and 
spectrometer  arrangement  especially  designed  for  Compton  scattering.  An  excellent  recent  review  of 
subsequent  measurements  and  theory  of  the  Compton  profile  has  been  given  by  Cooper  (1997). 

As  discussed  by  Duncanson  and  Coulson  (1945),  who  obtained  closed-form  analytical  expressions  for 
the  Compton  line-shapes  for  isolated  (unbound)  atoms  He  through  Ne,  if  a monoenergetic  photon 
beam  falls  on  an  atom  in  which  all  the  electrons  are  at  rest,  the  scattered  photons  at  any  given 
deflection  angle  should  have  the  single  energy  predicted  by  the  Compton  (1922,  1923)  and  Debye 
(1923)  formula.  However,  if  the  electrons  are  in  motion,  as  we  know  to  be  the  case,  there  is  a 
Doppler  effect  related  to  the  projected  velocities  of  the  electrons.  This  causes  the  modified 
(Compton)  line  to  spread  out  and  become  a band,  whose  profile  can  then  be  used  to  determine  the 
electron  velocity  distribution.  This  Doppler-based  description  (Du  Mond  1933)  gives  effectively  the 
same  results  as  the  complete  wave-mechanical  theory  of  Waller  (1927)  and  Wentzel  (1927a,b). 

The  newly  available  high  intensity  and  higher-energy  synchrotron  radiation  sources,  heading  toward 
the  50  to  500  keV  range,  according  to  Cooper  (1997),  make  possible  the  exploitation  of  Compton 
profile  measurements  to  examine  the  wave-mechanical  structures  of  the  scatterers.  For  example, 
Compton  profiles  can  be  used  to  check  band  theoretical  models  of  solids,  to  identify  electron 
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correlation  effects,  and  recently  have  been  used  to  extract  Fermi  surface  topography  in  Li  metal 
(Schiilke  et  al  1996). 

Clearly,  high-resolution  Compton  profile  measurements  are  a useful  tool  for  a variety  of  atomic 
structure  and  solid  state  investigations.  However,  the  impact  of  our  knowledge  of  these  profiles  on 
general-use  systematic  compilations  of  x-ray  attenuation  coefficients,  pip,  such  as,  for  example,  those 
of  Berger  and  Hubbell  (1987,  1998/1999),  Creagh  and  Hubbell  (1992),  Hubbell  and  Seltzer  (1995) 
•and  of  Cullen,  Hubbell  and  Kissel  (1997),  also  on  systematic  calculations  and  compilations  of  the 
incoherent  scattering  cross  section  such  as  in  Hubbell  et  al  (1975),  has  not  yet  been  specifically 
treated  in  the  literature,  to  this  author’s  knowledge.  It  is  the  aim  of  this  report,  in  the  survey  in 
Section  3,  and  in  the  conclusions  in  Section  4,  to  provide  a basis  for  such  a treatment. 


2.  Survey  of  S-Matrix  Theory. 

A concise  definition  of  the  S-matrix  (scattering  matrix,  also  sometimes  called  the  collision  matrix)  is 
given  in  The  International  Dictionary  of  Physics  and  Electronics  (Michels  1961):  "Consider  a 
quantum  mechanical  system  composed  of  two  parts,  such  as  (1)  an  atom  or  nucleus,  and  (2)  the 
electromagnetic  radiation  field,  or  such  as  (1)  an  incoming  proton  and  (2)  a nucleus.  It  is  possible  to 
express  the  initial  state  of  the  system  as  a linear  combination  of  eigenfunctions.  After  a photon  of  the 
electromagnetic  field  has  been  scattered  from  the  atom  or  nucleus,  or,  for  example  in  the  case  of 
particle-radiation  interactions,  after  a proton  has  been  scattered  from  a nucleus,  the  final  state  of  the 
system  can  again  be  expressed  as  a different  linear  combination  of  the  same  set  of  eigenfunctions  as 
that  used  in  the  initial  description.  If  the  wave  functions  for  the  initial  and  final  states  are  ^(t  = -oo) 
and  rp(t  = + oo ) respectively,  and  are  given  by  the  expressions 


= -os)  = EtAA,  and  (1) 

H t = +o°)  = Embm^m,  (2) 

in  terms  of  the  system’s  eigenfunctions  \[/B,  then  the  matrix  having  elements  S^,  such  that 
bm  = ^n^nuAi  (3) 

is  called  the  S-matrix. " 

According  to  Newton  (1982),  the  S-matrix  was  first  independently  introduced  in  quantum  scattering 
theory  by  Wheeler  (1937)  and  by  Heisenberg  (1943  a,b),  and  is  sometimes  referred  to  as  the  S matrix 
of  Heisenberg  (Goldberger  and  Watson  1964).  Treatments  of  the  S-matrix  and  its  applications  appear 
in  many  standard  textbooks,  for  example  Blatt  and  Weisskopf  (1952),  Ballentine  (1998)  and  in 
particular  Rodberg  and  Thaler  (1967)  which  devotes  an  entire  chapter  to  the  S-matrix  and  the  K- 
matrix,  the  latter  of  which  conveys  the  same  information  as  the  S-matrix  but  has  different  properties. 

The  first  significant  application  of  the  S-matrix  formalism  to  Compton  scattering  appears  to  be  that  by 
Whittingham  (1971,  1981)  in  which,  according  to  Suric  et  al.  (1991),  the  external  field  was  assumed 
to  be  point  Coulombic  (unscreened),  and  a limited  number  of  partial  waves  and  multipoles  were 
included  in  the  calculations.  Whittingham’ s results,  according  to  Suric  et  al.  (1991),  were 
significantly  larger  than  the  predictions  of  more  approximate  methods,  for  the  cases  that  were 
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considered,  and  the  computational  codes  could  not  be  verified.  Another  early  application  of  S-matrix 
formalism  to  Compton  scattering  was  that  by  Wittwer  (1972). 

As  yet,  there  appear  to  be  no  systematic  (all-Z,  all-energies)  computations  and  tabulations  of  Compton 
cross  sections,  where  Z is  the  atomic  number,  utilizing  the  S-matrix  formalism,  in  a form  amenable  to 
incorporation  in  systematic  (all-Z,  extended  photon  energy  range)  photon  cross  section  and  attenuation 
coefficient  data  bases.  The  existing  available  theoretical  tools  for  improved  Compton  scattering 
computations,  including  "exact"  second-order  relativistic  S-matrix  prescriptions,  have  recently  been 
reviewed  by  Bergstrom  and  Pratt  (1997).  This  review  outlines  the  recent  Compton  theory  progress  to 
date,  primarily  from  the  efforts  of  a series  of  collaborations  between  researchers  at  the  University  of 
Pittsburgh,  the  Lawrence  Livermore  National  Laboratory,  and  the  Ruder  Boskovic  Institute  in  Zagreb, 
Croatia  (Suric  et  al.  1991,  Suric  1992,  Bergstrom  et  al.  1992,  Bergstrom  1993a,  Bergstrom  et  al. 
1993b,  Pratt  et  al.  1994,  Bergstrom  1995  and  Suric  et  al.  1996).  In  this  series  of  collaborations,  an 
exact  second-order  S-matrix  code  has  been  developed  for  the  relativistic  numerical  calculation  of  cross 
sections  for  Compton  scattering  of  photons  by  bound  electrons  within  the  independent  particle 
approximation  (IPA).  In  the  sample  cases  studied,  comparison  with  measurements  is  found  to  be 
better  than  with  impulse  approximation  calculations. 

The  recent  review  by  Kane  (1997)  of  experimental  studies  of  inelastic  (Compton)  scattering  of  x-ray 
and  y-ray  scattering  by  atoms  includes  some  comparisons  with  S-matrix  computations  as  well  as  with 
the  non-relativistic  Brown/Cromer  incoherent  scattering  factor  S(x,Z)  systematic  compilation  of 
Hubbell  et  al.  (1975).  Because  the  quantity  S(x,Z)  was  introduced  in  the  x-ray  crystallography 
literature,  the  momentum  transfer  parameter  x is  customarily  expressed  in  terms  of  the  photon 
wavelength  X in  Angstroms  (A),  where  1 A = 10'10  m,  rather  than  in  terms  of  the  photon  energy  E 
(e.g.,  in  keV,  where  1 A = 12.3984244/E(keV),  Cohen  and  Taylor  1986,  1997)  for  tabulation 
purposes.  Thus  in  the  definition:  x = sin(0/2)/A(A),  where  8 is  the  scattered-photon  deflection  angle, 
x has  the  units  A'1.  In  the  example  in  Kane  (1997)  for  Sn  (Z  = 50),  the  measurement  error  bars  for 
a Compton  profile  (279.2  keV,  8 = 115°,  where  0 is  the  scattered-photon  deflection  angle)  are  too 
large  to  draw  firm  conclusions  toward  preferring  the  S-matrix  results  over  the  non-relativistic  widely- 
used  incoherent  scattering  function  values.  The  comparisons  by  Namito  et  al.  (1995)  of  recent 
measurements  of  S(x,Z)  with  the  Hubbell  et  al.  (1975)  compilation  suggest  that  the  latter  compilation, 
in  addition  to  its  convenience,  is  still  useful  for  most  practical  applications. 


3.  Survey  of  Selected  Compton  Profile  Measurements  and  Theoretical  Treatments. 

The  question  has  been  raised,  of  whether  incorporation  of  the  effect  of  Doppler  broadening  on  the 
spectrum  of  the  Compton-scattered  photon  should  have  any  significant  effect  on  the  calculation  of  the 
cross  section  for  this  process,  for  providing  more-accurate  general-purpose  systematic  tabulations  of 
the  mass  attenuation  coefficient  pip  and  related  quantities  such  as  the  mass  energy-absorption 
coefficient  pjp.  If  one  is  interested  in  the  spectral  distribution  of  the  Compton  scattered  photons, 
and  not  just  the  integrated  cross  section,  Namito  et  al.  (1994)  have  made  it  dramatically  clear,  by 
computing  the  scattered  spectrum  with  and  without  including  Doppler  broadening  in  the  EGS4  code 
(Nelson  et  al.  1985),  that  such  broadening  must  be  included  in  order  to  replicate  the  corresponding 
measurements  by  Namito  et  al.  (1994)  using  40  keV  polarized  photons  from  a 2.5  GeV  synchrotron 
light  source  (KEK-PF). 

This  section  presents  a more-or-less  chronological  survey  and  account  of  selected  experimental  and 
theoretical  studies  of  the  Compton  profile,  or  Doppler  spectral  broadening  of  what  would  be  a 
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monoenergetic  scattered  photon,  for  a monoenergetic  source  and  a particular  photon  deflection  angle 
9.  Whether  or  not  any  conclusions  can  be  drawn  from  this  account,  at  the  very  least  this  section  (and 
the  references  cited)  provides  a guide  to  on-file  documents  on  this  topic  that  are  held  in  the  NIST 
Photon  and  Charged  Particle  Data  Center  reference  collection. 

The  first  major  theoretical  treatment  of  the  energy-broadening  of  the  Compton  scattered  line  appears 
to  be  that  of  Jauncey  (1925a,b,  1934)  who  considered  the  kinetic  energy,  and  hence  its  component 
distribution  of  velocities,  of  the  electron  in  its  Bohr  orbit.  The  next  step  was  the  work  by  Waller 
(1927)  who  used  the  Kramers-Heisenberg  dispersion  formula  to  show  that  one  should  see  a slightly 
broadened  Compton  line.  Throughout  the  1930s,  extensive  experimental  and  theoretical  studies  of  the 
Compton  line  structure  were  carried  out  by  Du  Mond  and  Kirkpatrick  (Du  Mond  1929,  1930,  Du 
Mond  and  Kirkpatrick  1931,  1937,  Kirkpatrick  and  Du  Mond  1938);  much  of  this  work  was 
organized  and  presented  in  an  extensive  and  definitive  Reviews  of  Modem  Physics  article  (Du  Mond 
1933).  As  with  Jauncey’s  earlier  work,  Du  Mond’s  careful  measurements  could  to  a good 
approximation  be  interpreted  in  terms  of  Doppler  broadening  due  to  the  plus  and  minus  velocity 
projections  of  the  orbiting  electrons. 

During  this  same  period,  Kappeler  (1936)  made  Compton  profile  measurements  on  gas  targets  (Ne, 

02  and  N2)  and  on  solids  (Li,  C and  NaF),  and  Burkhardt  made  measurements  on  Ne.  Schnaidt 
(1934)  and  Franz  (1937)  added  refinements  to  the  theoretical  understanding  of  the  process.  Hicks 
(1937)  used  self-consistent-field  variation  functions  which  take  into  consideration  the  instantaneous 
interaction  of  the  electrons,  momentum  distributions  and  intensity  distributions  to  compute  Compton 
lines  for  helium  and  molecular  hydrogen.  Thes  computations  neglected  the  small  relativity  and 
binding  corrections. 

Moving  into  the  World  War  II  period,  Duncanson  and  Coulson  (1945)  derived  closed  analytical 
expressions  for  the  shapes  of  the  Compton  lines  scattered  from  the  elements  H to  Ne,  in  the  form  of 
free  atoms.  From  their  analysis  of  the  then-available  body  of  measurements,  they  concluded  that  by 
forming  a diatomic  molecule,  the  width  of  the  Compton  profile  is  increased  by  about  15%  to  25%. 

Brini  et  al.  (1960)  measured  the  differential  Compton  cross  section  for  662  keV  photons  on  the  K- 
shell  electrons  of  Pb,  using  the  coincidence  technique  with  the  K x rays  produced  in  the  process,  and 
with  care  to  use  Cu  shielding  to  avoid  detecting  K x rays  from  the  Pb  collimation  system.  In  this 
same  period,  Kilby  experimentally  investigated  the  electron  momentum  distributions  in  H2  and  He 
(Kilby  1961)  and  in  metallic  Li  (Kilby  1963),  and  concluded  that  there  were  some  serious 
discrepancies  between  theory  and  experiment  with  respect  to  the  mean  kinetic  energy  of  the 
conduction  electrons,  at  least  in  Li.  Motz  and  Missoni  (1961)  and  Di  Lazzaro  and  Missoni  (1966) 
also  measured  Compton  spectra  from  scattering  by  K-shell  electrons  of  Sn,  Au  and  U,  and  Varma  and 
Eswaran  (1962)  measured  Compton  spectra  from  the  K-shell  electrons  of  Pb.  All  measurements  were 
made  with  662  keV  photons. 

Kilby  (1965)  derived  expressions  for  the  Compton  profile  intensity  based  on  the  wave-mechanical 
scattering  theory  of  Waller  and  Hartree  (1929)  assuming:  (a)  a one-electron  approximation  is  valid, 

(b)  electronic  binding  energies  are  small  compared  with  the  energy  transferred  from  the  photon,  and 

(c)  the  electron  wave  functions  for  the  relevant  continuum  states  are  approximated  by  plane  waves. 

He  concluded  from  this  study  that  the  shape  of  the  Compton  profile  near  its  peak  is  particularly 
sensitive  to  the  description  of  the  more  loosely  bound  electrons,  and  that  hence  the  Compton  profile 
measurements  can  yield  information  about  the  valence  electrons  of  the  heavier  elements.  Veigele  et 
al.  (1966,  1969)  used  a relativistically  covariant  formalism  to  derive  a generalized  Compton  equation; 
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the  derived  Compton  profile  showed  a defect  (a  reduced  intensity)  in  the  broadened  line  toward  a 
wavelength  shorter  than  that  predicted  by  the  Compton  equation. 

Weiss  and  Phillips  and  collaborators  (Weiss  et  al.  1968a,  Weiss  and  Phillips  1968b,  Phillips  and 
Weiss  1968,  Weiss  1972,  1973a,b,  1975,  Phillips  and  Weiss  1972a,b,  Phillips  1973,  Weiss  et  al. 
1977)  exploited  Compton  profiles  for  examining  valence  electrons  in  solids,  both  experimentally  and 
theoretically,  and  published  extensive  tables  of  Compton  line  shapes  and  electron  momentum 
distributions  for  atoms  and  ions  as  heavy  as  Ge  (Z=32).  Schulke  et  al.  (1969)  measured  directional 
Compton  profiles  for  Be,  and  found  evidence  for  Coulomb  interactions  between  electrons  that  were 
consistent  with  the  RPA  (random  phase  approximation)  calculations  of  Ohmura  and  Matsudaira  (1964) 
for  photon  inelastic  scattering  by  an  electron  in  a degenerate  electron  gas.  In  turn,  Ohmura  and 
Matsudaira  (1964)  noted  a discrepancy  between  the  Kappeler  (1936)  Li  Compton  profile  measured 
width  and  the  Wigner-Seitz  wave  function  calculation  by  March  (1954)  and  Donovan  and  March 
(1956)  which  could  not  be  explained  by  Coulomb  correlation.  More  recently,  Schulke  et  al.  (1996) 
used  synchrotron  radiation  from  the  DORIS  storage  ring,  monochromatized  to  31  keV,  to  measure 
directional  Compton  profiles  from  Li  metal.  They  found  discrepancies  from  the  recent  density- 
functional  calculations  by  Sakurai  et  al.  (1995)  that  they  attributed  to  self-energy  effects  connected 
with  the  excitation  of  so-called  plasmaron  modes. 

Another  extensive  body  of  original  and  review  work  on  Compton  profiles  is  that  by  Cooper  and 
collaborators  (Cooper  et  al.  1970,  1974,  1976,  Stuewer  and  Cooper  1977,  Cooper  1985,  Cardwell 
and  Cooper  1986,  Cooper  1997).  In  one  of  these  works,  Cooper  (1997)  points  out  that  what  has  been 
called  in  common  parlance  "magnetic  Compton  scattering"  should  now  be  called  "spin-dependent 
Compton  scattering."  During  this  same  period,  Eisenberger  and  collaborators  (Eisenberger  1970, 
1972a,  Eisenberger  et  al.  1972b, c,d,  Reed  and  Eiesenberger  1972,  1974)  also  have  provided 
extensive  Compton  profile  measurements  and  comparisons  with  theory.  Their  observations  of  high- 
momentum  anisotropy  in  both  Cu  and  Ni  suggested  d-electron  bonding  in  the  systems  studied. 
Manninen  and  collaborators  (Paakkari  et  al.  1972,  1973a,b,  Manninen  1974,  Paatero  et  al.  1974, 
Manninen  1978)  developed  an  annular  collimation  arrangement  with  which  they  made  Compton 
profile  measurements  with  Al,  V and  Fe  targets.  They  found  considerable  deviations  in  their  results 
from  predictions  of  the  free-atom,  free-electron  model  at  low  momentum  transfers.  However,  there 
was  reasonable  agreement  for  high  momentum  transfers. 

Epstein  (1973)  measured  Compton  profiles  in  benzene  and  several  other  hydrocarbons,  and  Tanner 
and  Epstein  (1974)  continued  these  studies  with  Compton  profiles  from  a water  target.  From  a 
comparison  of  their  measurements  with  "near  Hartree-Fock"  (NHF)  and  "minimal  basis  set"  (MBS) 
calculations,  they  concluded  that  neither  energy  minimization  nor  fulfillment  of  the  virial  theorem  is 
in  itself  sufficient  to  guarantee  a reliable  momentum  distribution. 

Ribberfors  (1975a)  developed  expressions  for  Compton  profiles  in  the  case  of  isotropic  momentum 
distributions,  and  generalized  these  results  (Ribberfors  1975b)  to  include  the  effects  of  both  anisotropy 
and  polarization.  In  a further  work,  Ribberfors  and  Berggren  (1982)  showed  that  incoherent 
scattering  functions  S(x,Z),  defined  as  an  integral  over  Compton  profiles,  could  easily  ("on  a pocket 
calculator")  be  computed  by  means  of  a linear  approximation.  However,  their  results  for  low  Z 
materials  were  barely  distinguishable  from  the  non-relativistic  S(x,Z)  compilation  of  Hubbell  et  al. 
(1975),  both  being  well  within  the  claimed  uncertainties  of  available  measurements,  and  for  a high  Z 
material  (e.g.,  Pb)  their  results  were  completely  indistinguishable.  On  the  other  hand,  in  a recent 
personal  communication,  Ribberfors  (1999)  cautioned  that  for  low  Z and  low  x there  could  be  a 
substantial  difference.  For  example,  Ribberfors  and  Berggren  (1982),  found  for  Al  and  x < 0.2  A'1 
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that  the  integrated  Compton  profile  S(x,Z)  results  were  as  much  as  40%  less  than  the  Waller-Hartree 
non-relativistic  S(x,Z)  Cromer  and  Mann  (1967b)  and  Cromer  (1969)  values  incorporated  into  the 
Hubbell  et  al.  (1975)  compilation. 

Williams  et  al.  (1976),  acting  in  behalf  of  the  International  Union  of  Crystallography  Commission  on 
Spin,  Charge  and  Momentum  Density,  endeavored  to  standardize  Compton  profile  measurements, 
taking  water  as  their  anchor  material  for  this  project.  Five  different  laboratories  participated, 
providing  13  independent  Compton  profile  measurements  on  water.  Inconsistent  normalization  of  the 
various  measured  profiles  was  found  to  be  the  major  source  of  error.  After  renormalization  of  the 
profiles,  the  results  were  found  to  be  consistent  to  within  ±2  % of  the  peak  height  of  the  profiles.  It 
was  concluded  that  reproducibility  of  ± 1 % of  the  peak  height  of  the  profile  is  quite  feasible,  and 
that  in  a differential  experiment  ±0.5  % or  even  better  is  possible.  The  techniques  are  extensively 
reviewed  in  Williams  (1977b),  with  further  details  in  the  chapters  in  Williams  (1977a)  by  Halonen  et 
al.  (1977),  Mendelsohn  and  Smith  (1977)  and  by  Berggren  et  al.  (1977). 

Owen  (1977,  1978)  theoretically  examined  both  binding  effects  and  "off-mass-shell"  effects  on  the 
Compton  profile  using  Dirac  wave  functions.  From  his  calculations,  he  concluded  that  ignoring  these 
effects  for  photons  scattered  in  the  forward  directions  from  electrons  in  high-Z  atoms  distorts  and 
extends  the  low-energy  part  of  the  profile.  Further,  he  observed  that  these  uncorrected  profiles  have 
secondary  maxima  which  contribute  a large  spurious  amount  to  the  calculated  cross  section. 

Pattison  and  Schneider  (1979)  made  a test  of  the  relativistic  Is  wave  functions  in  Au  and  Pb  by 
measuring  the  Compton  profiles  in  these  elements  using  412  keV  gamma  radiation.  Their  observed 
energy  dependence  of  the  Compton  scattering  cross  section  was  in  good  agreement  with  the 
expression  given  by  Ribberfors  (1975a),  and  they  found  good  agreement  between  their  measurements 
and  the  relativistic  Hartree-Fock  Compton  profiles  calculated  by  Biggs  et  al.  (1975).  In  an  extension 
of  this  work,  Weyrich  et  al.  (1979)  examined  the  Compton  profiles  of  atoms  and  molecules  using 
Fourier  analysis.  In  particular,  they  investigated  the  possibility  of  extracting  accurate  bond  lengths, 
with  the  fluorine  molecule  F2  as  an  example. 

Carlsson  et  al.  (1982)  have  explored  the  inclusion  of  Compton  profile  effects  on  scattering  cross 
section  calculations  for  applications  in  diagnostic  medical  radiology.  They  provide  a useful 
assortment  of  formulas  and  procedures  for  this  purpose. 

Schumacher’s  group  (Reineking  et  al.  1983,  Wenskus  et  al.  1985)  calculated  double  differential  bound 
Compton  scattering  cross  sections  in  the  relativistic  impulse  aproximation  using  DFHS  (Dirac- 
Hartree-Fock-Slater)  wave  functions.  They  found  good  agreement  with  their  measurements 
(Schumacher  1971,  Rullhusen  and  Schumacher  1976,  Reineking  et  al.  1983)  of  Compton  scattering  of 
279  and  662  keV  photons  by  Cu,  Sn  and  Pb  targets  at  the  peak  of  the  profiles.  However,  on  the 
high-energy  sides  of  their  measured  profiles  they  noted  discrepancies  which  were  only  partly 
explained  by  second-order  S-matrix  calculations.  Kane  et  al.  (1987)  included  Compton  scattering 
measurements  of  88.03  keV  gamma  rays  on  a Bi  target  in  their  work  which  also  included  cross 
sections  for  elastic  scattering  and  atomic  photoeffect.  They  obtained  an  incoherent  scattering  value 
S(x,Z)  at  x = 6.3  A*1  of  68.6±4.7  to  compare  with  the  non-relativistic  value  of  77  in  the  Hubbell  et 
al.  (1975)  compilation.  This  work  is  included  in  the  extensive  review  by  Kane  (1992)  of  inelastic 
scattering  of  photons  in  the  x-ray  and  gamma-ray  energy  region. 

Most  recently,  Deb  et  al.  (1999)  have  measured  directional  Compton  profiles  for  indium  phosphide 
(InP)  using  gamma  rays  from  241  Am  and  compared  their  measured  results  with  their  calculated  results 
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which  were  based  on  the  linear  combination  of  Gaussian  orbitals  (LCGO)  method.  They  found  good 
agreement  and  that  the  anisotropy  dependences  have  an  intimate  connection  with  the  bonding 
properties  of  the  semiconductor.  This  work  was  a continuation  of  the  measurements  by 
Bandy opadhy ay  et  al.  (1996)  of  Compton  profiles  of  vanadium  carbide  and  theoretical  interpretation 
by  Deb  and  Chatteijee  (1996a)  and  measurements  by  Bandyopadhyay  et  al.  (1994)  of  profiles  of 
cuprous  oxide  and  theoretical  interpretation  by  Deb  and  Chatterjee  (1996b)  who  used  the  LCGO 
model  method.  Further  on  work  Compton  profile  measurements  and  theory  by  this  group,  including 
related  properties  such  as  density  of  states  (DOS)  and  charge  analysis,  with  emphasis  on  chemical 
bonding,  is  in  progress. 


4.  Conclusions  and  Recommendations. 

S-Matrix  Incoherent  Scattering  Theory: 

From  the  papers  examined  and  described  in  Section  2,  one  can  conclude  that  for  Compton  scattering 
of  photons  from  electrons  bound  in  atoms,  differential  and  integrated  cross  sections  for  this  process 
calculated  using  the  S-matrix  formalism  should  give  values  in  closer  agreement  with  measured  values. 
However,  there  is  a problem  in  presentation  of  the  results  for  practical  application,  such  as  in  EGS4 
codes  for  transport  calculations,  and  in  other  applications,  since  the  data  array  for  the  equivalent  of 
the  now  two-dimensional  incoherent  scattering  function  S(x,Z)  would  become  a three-dimensional 
function  S(E,0,Z)  in  which  E is  the  photon  incident  energy  and  d is  the  scattered-photon  deflection 
angle.  With  modem  computers,  such  higher  dimensional  arrays  are  perhaps  feasible. 

However,  as  noted  from  the  Namito  et  al.  (1995)  comparisons  of  incoherent  scattering  function 
measurements  with  theory,  the  non-relativistic  but  convenient  two-dimensional  S(x,Z)  values  compiled 
by  Hubbell  et  al.  (1975)  by  combining  and  extrapolating  the  works  of  Brown  (1970a,b,  1972,  1974), 
of  Cromer  and  Mann  (1967a,b)  and  of  Cromer  (1969,  1971)  appear  to  be  indistinguishable  from  the 
S-matrix  formalism  results,  within  the  error  bars  of  the  available  experiments,  particularly  for  high-Z 
elements.  Ribberfors  (1999)  points  out  the  significant  differences  between  non-relativistic  and 
relativistic  impulse-approximation  theoretical  treatments  for  Al  below  x = 0.2  A'1  (see.  e.g.,  Figure 
2 in  Ribberfors  and  Berggren  (1982)). 

Recently,  Kahane  (1998)  has  used  the  Ribberfors-Bergren  (1982)  treatment  to  compute  and  tabulate 
relativistic  Dirac-Hartree-Fock  photon  incoherent  scattering  functions  S(x,Z)/Z  for  all  elements  1 < 

Z < 110  for  a range  of  values  of  the  momentum  transfer  parameter  0.005  A'1  < x < 1000  A'1.  The 
only  measured  S(x,Z)  values  in  this  region,  with  which  to  compare  and  assess  non-relativistic  vs. 
relativistic  treatments,  are  those  of  Laval  (1942)  [in  France,  at  the  height  of  the  1940-1944  World 
War  II  Nazi  occupation]  for  which  no  error  estimates  are  provided;  also,  little  information  is  given  on 
the  experimental  details.  However,  in  the  future,  it  is  anticipated  that  relativistic  S(x,Z)  values  such 
as  those  by  Kahane  (1998)  will  replace  the  Hubbell  et  al.  (1975)  non-relativistic  S(x,Z)  values,  in 
both  incoherent  (Compton)  scattering  and  triplet  production  cross  section  computations  for  p/p 
compilation  purposes. 

Measured  and  Theoretical  Compton  Profile  Data: 

From  the  above  survey  in  Section  3 of  measured  and  theoretical  Compton  profile  data,  by  its  very 
nature  of  being  sensitive  to  chemical  bonding  by  valence  electrons,  it  would  be  very  difficult  to 
systematically  compile  "engineering  type"  tables,  such  as  is  now  done  for  the  mass  attenuation 
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coefficients  nip  and  incoherent  scattering  functions  S(x,Z),  for  scattering  and  attenuation  data  which 
depends  on  the  bonding-sensitive  Compton  profiles.  However,  as  dramatically  illustrated  by  Namito 
et  al.  (1994),  Doppler  broadening  should  and  can  be  incorporated  into  the  EGS4  (Nelson  et  al.  1985) 
and  other  Monte  Carlo  photon-electron  radiation  transport  codes,  as  shown  by  their  comparisons  of 
calculated  and  measured  spectra  for  C and  in  Cu. 
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